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The theoretical resuhs regarding the "transition frequencies" of two acoustically interact- 
ing bubbles have been verified numerically. The theory provided by Ida [Phys. Lett. A 297 
(2002) 210] predicted the existence of three transition frequencies per bubble, each of which 
has the phase difference of it/2 between a bubble's pulsation and the external sound field, 
while previous theories predicted only two natural frequencies which cause such phase shifts. 
Namely, two of the three transition frequencies correspond to the natural frequencies, while 
the remaining does not. In a subsequent paper [M. Ida, Phys. Rev. E 67 (2003) 056617], it 
was shown theoretically that transition frequencies other than the natural frequencies may 
cause the sign reversal of the secondary Bjerknes force acting between pulsating bubbles. 
In the present study, we employ a direct numerical simulation technique that uses the com- 
pressible Navier-Stokes equations with a surface-tension term as the governing equations to 
investigate the transition frequencies of two coupled bubbles by observing their pulsation am- 
plitudes and directions of translational motion, both of which change as the driving frequency 
changes. The numerical results reproduce the recent theoretical predictions, validating the 
existence of the transition frequencies not corresponding to the natural frequency. 

KEYWORDS: bubble dynamics, secondary Bjerknes force, direct numerical simulation, natural 
frequency, transition frequency 



1. Introduction 

The secondary Bjerknes force is an interaction force acting between pulsating gas bubbles 
in an acoustic field. The classical theory originated by Bjerknes predicts either attraction 
only or repulsion only, depending on whether the driving frequency stays outside or inside, 
respectively, the frequency region between the partial (or monopole) natural frequencies of 
two bubbles. However, recent studies show that the force sometimes reverses its own direction 
as the distance between the bubbles changes. The first theoretical study on this subject 
was performed by Zabolotskaya.^ Employing a linear coupled oscillator model, she showed 
that the radiative interaction between the bubbles, which results in the change in the natural 
frequencies of the bubbles, could cause this reversal. In the mid-1990s, Doinikov and Zavtrak 
arrived at the same conclusion by employing a linear theoretical model in which the multi- 
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pie scattering between the bubbles is more rigorously taken into account.^ These theoretical 
results arc considered to explain the stable structure formation of bubbles in a sound field, 
called "bubble cluster" or "bubble grape," which has been observed experimentally by several 
researchers in different fields.^' ''^^^ 

In both of the theoretical studies mentioned above, it was assumed that the reversal is 
due to the change in the natural (or the resonance) frequencies of bubbles, caused by the 
radiative interaction between bubbles. However, those authors had differing interpretations 
of how the natural frequencies change. The theoretical formula for the natural frequencies, 
used by Zabolotskaya'^ and given previously by Shima,^^ shows that the higher and the lower 
natural frequencies (converging to the partial natural frequencies of a smaller and a larger 
bubble, respectively, when the distance between the bubbles is infinite) reveal an upward and 
a downward shift, respectively, as the bubbles come closer to one another. In contrast, 
Doinikov and Zavtrak assumed intuitively that both the natural frequencies rise.^'^ This 
assumption seems to explain well the sign reversal occurring not only when both bubbles are 
larger than the resonance size but also when one bubble is larger and the other is smaller 
than the resonance size. The sign reversal in the latter case, for instance, is thought to occur 
when the resonance frequency of a larger bubble, increasing as the bubbles come closer to 
one another, surpasses the driving frequency, resulting in the change in the pulsation phase 
of the larger bubble, and leading to the change in the phase difference between the bubbles. 
However, this assumption is obviously inconsistent with the previous theoretical result for the 
natural frequencies.^' 

Recently, Ida^'^ proposed an alternative theoretical explanation for this phenomenon, also 
using the linear model Zabolotskaya used. He claimed that this phenomenon cannot be inter- 
preted by only observing the natural frequencies, and that it is useful to define the transition 
frequencies that make the phase difference between a bubble's pulsation and an external sound 
7r/2 (or 37r/2). It has been pointed out theoretically that the maximum number of natural 
frequencies and that of transition frequencies are, in general, not in agreement in multibubble 
cases, ^^'^^ while they are, as is well known, consistent in single-bubble cases where the phase 
difference between a bubble's pulsation and an external sound becomes 7r/2 only at its natural 
frequency. (This is not true in strongly nonlinear cases, in which the phase reversal can take 
place even in frequency regions far from the bubble's natural frequency; see, e.g., refs. 16,17.) 
In a double-bubble case, for instance, that theory predicts three transition frequencies per 
bubble, two of which correspond to the natural frequencies.^^ A preliminary discussion on a 
A?^-bubble system^^ showed that a bubble in this system has up to 2N — 1 transition frequen- 
cies, only N ones of which correspond to the natural frequencies. More specifically, the number 
of transition frequencies is in general larger than that of natural frequencies. The transition 
frequencies not corresponding to the natural frequency have differing physical meanings from 
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those of the natural frequencies; they do not cause the resonance response of the bubbles. "^^'"^^ 
The theory for the sign reversal of the force, constructed based on the transition frequencies, 
predicts that the sign reversal takes place around those frequencies, not around the natural 
frequencies, and can explain the sign reversal in both cases mentioned above. -"^^ Moreover, the 
theory does not contradict the theory for the natural frequencies described previously, because 
all the natural frequencies are included in the transition frequencies. 

The aim of this paper is to verify the theoretical prediction of the transition frequencies by 
direct numerical simulation (DNS). In a recent paper, Ida proposed a DNS technique, based 
on a hybrid advection scheme, ""^^ a multi-time-step integration technique, ""^^ and the Cubic- 
Interpolated Propagation/Combined, Unified Procedure (CIP-CUP) algorithm, ^'^'^^ which 
technique allows us to compute the dynamics (pulsation and translational motion) of de- 
formable bubbles in a viscous liquid even when the separation distance between the bubbles is 
small. -"^^'^^ In that DNS technique, the compressible Navier-Stokes equations with a surface- 
tension term are selected as the governing equations, and the convection, the acoustic, and 
the surface-tension terms, respectively, in these equations are solved by an explicit advection 
scheme employing both interpolation and extrapolation functions which realizes a discontin- 
uous description of interfaces between different materials, ""^^ by the Combined, Unified Pro- 
cedure (CUP) being an implicit finite difference method for all-Mach-numbcr fiows,^'^'^^ and 
by the Continuum Surface Force (CSF) model being a finite difference solver for the surface- 
tension term as a volume force. Efficient and accurate time integration of the compressible 
Navier Stokes equations under a low-Mach-mmibcr condition is achieved by the multi-timc- 
step integration technique, which solves the different-nature terms in these equations with 
different time steps. Further details of this DNS technique can be found in refs. 18, 19. Em- 
ploying this DNS technique, in the present study we perform numerical experiments involving 
two acoustically coupled bubbles in order to investigate the recent theories by observing the 
bubbles' pulsation amplitudes and the directions of their translational motion. In particular, 
we focus our attention on the existence of the transition frequencies that do not cause the 
resonance response. 

The following sections are organized as follows. In §2, the previously expounded theories 
are reviewed and reexamined, including those for the transition frequencies and the sign 
reversal, and in §3, the numerical results and a discussion are provided. Section 4 presents 
concluding remarks. 

2. Theories 

2.1 A single-bubble problem 

It is well known that, when the wavelength of an external sound is sufficiently large com- 
pared to the radius of a bubble (named "bubble 1," immersed in a liquid) and the sphericity 
of the bubble is maintained, the following second-order differential equation^' describes the 
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linear pulsation of the bubble: 

ei + WiqCi + diei = — , (1) 

Po^io 



1 



3kPo + (3ac-1)|^ 



where it is assumed that the bubble's time-dependent radius can be represented by i?i = 
Rio + ei (|ei| ^ -Rio) with Riq being its equilibrium radius and ei being the deviation in the 
radius, and ljiq is the bubble's natural frequency, is the damping factor, ^^'^^ pcx is the sound 
pressure at the bubble position, po is the equilibrium density of the liquid, Pq is the static 
pressure, k is the polytropic exponent of the gas inside the bubble, a is the surface tension, 
and the over dots denote the derivation with respect to time. Assuming pcx = —PaSvuuot 
(where Pa is a positive constant), the harmonic steady-state solution of eq. (1) is determined 
as 

ei = Ksism.{ut - (j)si), 

where 

Pa I 1 



Ksi = 



1 f Siu; 



(psi = tan" 



UJ 



This result reveals that the phase reversal (or the phase difference of (j)si = tt/2) appears at 
only to = uJiQ. Moreover, if 6i <^ uio, the bubble's resonance response can take place at almost 
the same frequency. Though the resonance frequency shifts away from ujiq as 6i increases, it 
can in many cases be assumed to be almost the same as cuiq. (Also, the nonlinearity in bubble 
pulsation is known to alter a bubble's resonance frequency.^' ^'^'') 



2.2 A double-bubble problem 

When one other bubble (bubble 2) exists, the pulsation of the previous bubble is driven 
by not only the external sound but also the sound wave that bubble 2 radiates. Assuming 
that the surrounding liquid is incompressible and the time-dependent radius of bubble 2 can 
be represented by R2 = R20 + ^2 (|e2| <C -R20)) the radiated pressure field is found to be,^^ 

p{r,t) ^ -62, 

r 

where r is the distance measured from the center of bubble 2. The total sound pressure at the 
position of bubble 1 (pdi) is, thus, 

P^l^Pe. + ^e2, (2) 
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where D is the distance between the centers of the bubbles. This total pressure drives the 
pulsation of bubble 1. 

Replacing in eq. (1) with p^i yields the modified equation for bubble 1, 

•• I 2 I X • ^'ex -^20 •• /o^ 

ei+a;ioei + ()iei = -— TT^^a, (3) 

and exchanging 1 and 2 (or 10 and 20) in the subscripts in this equation yields that for bubble 
2, 

62 + i02Qe2 + 6262 = -— ^-^ei. (4) 

This kind of system of differential equations is called a (linear) coupled oscillator model or 
a self-consistent model,^^'^^ and is known to be third-order accuracy with respect to 
The same (or essentially the same) formula has been employed in several studies considering 
acoustic properties of coupled bubbles.^' ^^'^^"^^ 

Shima-*^^ and Zabolotskaya,^ assuming Sj pa (for j = I and 2), derived the theoretical 
formula for the natural frequencies, represented as 

/ 2 2\/ 2 2\ -R1O-R2O 4 n 

(cJiO )(W20 ) ^2— W 

This equation predicts the existence of up to two natural frequencies (or eigenvalues of the 
system (3) and (4)) per bubble. Exchanging 10 and 20 in the subscripts in this equation 
yields the same equation; namely, when 5j ~ 0, both the bubbles have the same natural 
frequencies. The higher and the lower natural frequencies reveal an upward and a downward 
shift, respectively, as D decreases. ""^^'^^ 

The theoretical formula for the transition frequencies, derived by Ida,^'^ can be obtained 
based on the harmonic steady-state solution of system (3) and (4). Assuming pex = —Pa sin cut, 
the solution is determined as follows: 



ei = Ki sin(a;t — 0i), (6) 



where 



with 



01 = tan-i (8) 



_ HiF + M2G _ HiG - M2F 

F = L1L2 ^2 — ^ ~ ^1^2, 

G = L1M2 + L2M1, Hi = L2 + 
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(b) Transition frequencies for S.~0 
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Fig. 1. Transition frequencies as functions of the distance between the bubbles. The lower figure 
shows the transition frequencies in the case where the damping effect is neglected. The dashed 
lines denote the transition frequencies that do not cause resonance. The arrows are used for a 
discussion described in the text. 



-^1 — '^10 ~ ^ ) L2 — ^20 ~ ; 
Ml = SlLO, M2 = d2UJ. 

Exchanging 1 and 2 (or 10 and 20) in the subscripts in these equations yields the formula for 
bubble 2. Based on the definition, the transition frequencies are given by^^ 

HiF + M2G = 0. (9) 

This equation predicts the existence of up to three transition frequencies per bubble; this 
number is greater than that of the natural frequencies given by eq. (5) . This result means that 
in a multibubble case, the phase shift can take place not only around the bubbles' natural 
frequencies but also around some other frequencies. Moreover, it was shown in ref. 13 that the 
transition frequencies other than the natural frequencies do not cause the resonance response; 
namely, one of the three transition frequencies has different physical meanings from those of 
the natural frequency. 

Figure 1(a) shows the transition frequencies of bubbles 1 and 2, uji and UJ2, for Rio = 5 
fim and R20 = 9 /xm, as functions of D/{Rio + i?2o)- These bubble radii are chosen to be small 
enough so that the sphericity of the bubbles is maintained sufficiently. The parameters are 
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set to pQ = 1000 kg/m^, k = 1.33, Pq = 1 atm, and a = 0.0728 N/m. The damping factors 
are determined by the sum of the viscous and radiation ones as Sj = (4/i/po^jo) + ("^^-^jo/c) 
since the DNS technique does not consider the thermal conduction/^ where the viscosity of 
the hquid fi = 1.137 x 10~^ kg/(m s) and its sound speed c = 1500 m/s. As a reference, 
in Fig. 1(b) we display the transition frequencies for 6j 0. As discussed previously, the 
highest and the second highest transition frequencies of the larger bubble tend to vanish when 
the damping effect is sufficiently strong; in the present case, they disappear completely. The 
second highest and the lowest transition frequencies of the smaller bubble cross and vanish 
at a certain distance, and only one transition frequency remains for sufficiently large D. For 
D ^ oo, the transition frequencies remaining converge to the partial natural frequencies of the 
corresponding bubble. The solid lines displayed in Fig. 1(b) denote the transition frequencies 
that correspond to the natural frequencies given by eq. (5). 

We clarify here physical meanings of the transition frequencies that do not accompany 
resonance, which meanings have not yet been described in the literature. Let us consider 
the bubbles under the condition indicated by the dots (the origin of the arrows) shown in 
Fig. 1. Bubble 2 under this condition emits a strong sound (denoted below by p2{D)) whose 
oscillation phase is (almost) out-of-phase with the external sound pcxj because the driving 
frequency stays near and slightly above the natural frequency of this bubble. If P2{D) is 
measured at a point sufficiently near bubble 2, its amplitude can be larger than that of pex 
and hence the phase of the total sound pressure, pcx +P2{D), may be almost the same as that 
of p2{D), i.e., almost out-of-phase with pcx- As the driving frequency shifts toward a higher 
range along the arrows, the absolute value of P2{D) decreases and, at a certain frequency, 
becomes lower than |pcx|; the phase of the total sound pressure finally becomes (almost) in- 
phase with pcx- This transition of the power balance of the two sounds results in the phase 
reversal of bubble 1 without accompanying the resonance response. 

2. 3 The secondary Bjerknes force 

The secondary Bjerknes force acting between two pulsating bubbles is represented by^"^ 

Foc{ViV2) ,'"'~\ , (10) 
||r2 - rill"" 

where Vj and denote the volume and the position vector, respectively, of bubble j, and 
(■ ■ ■ ) denotes the time average. Using eqs. (6)-(8), this equation can be rewritten as^ 

F oc KiK2 cos(0i - </'2)-^^— (11) 

||r2 - rill 

The sign reversal of this force occurs only when the sign of cos(^i — ^2) (or of {V1V2)) changes, 
because i^i > and K2 > 0. Namely, the phase property of the bubbles plays an important 
role in the determination of the sign. Roughly speaking, the force is attractive when the 
bubbles pulsate in-phase with each other, while it is repulsive when they pulsate out-of-phase. 
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Fig. 2. Sign of the secondary Bjerknes force [cos{(j)\ — ^2)] and phase differences and (^2], deter- 
mined theoretically, as functions of the driving angular frequency. The positive value of cos((j!)i — ^2) 
indicates the attraction, while the negative one indicates the repulsion. The dashed lines denote 
the results for 5j w 0. 

The solid lines displayed in Fig. 2 denote cos{(f)\ — <j)2)-, (t^i, and (f)2 as functions of uj/loio, 
and the dashed lines denote those for Sj 0. The physical parameters are the same as used 
for Fig. 1, except for the separation distance fixed to I? = 20 /xm [D/{Rio + R20) ~ 1.43]. In 
this figure, we can observe sign reversals of cos(0i — ^2) at u/uio pa 0.72 and uj/loiq ~ 1.15. 
In ref. 14, it was shown theoretically that both the sign reversals are due to the transition 
frequencies that do not correspond to the natural frequencies; that is, the sign reversals take 
place near (or, when Sj — > 0, at) those frequencies. The present result follows that theoretical 
prediction. The respective reversals are observed near the second highest transition frequency 
of bubble 1 and near the highest of bubble 2 for Sj — > 0. Meanwhile, we can observe <^2 > tt, 
which should not be observed in a single-bubble case. Ida^^ explained that such a large phase 
delay is realizable in a multibubble case, by the radiative interaction. 

We make here a remark regarding the estimation of the points at which the sign reversal 
of the force takes place. It was shown previously^^ that the transition points of the force 
are hardly changed by the damping effects, even when bubbles are small {Rjo ~ 1 /xm) or 
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relatively large {Rjo ^ 1 mm). The present result displayed in Fig. 2 follows that finding. 
These results may allow us to consider that the simple theoretical formulas for the transition 
frequencies not causing resonance, 

u;l = ^ ^ and = ^ ^ 

^ I-R20/D 2 l-Rio/D 

derived for (5i — > and 62 O/^ can be a good approximation of the transition points, even 

in the cases of Sj 7^ 0. These formulas yield (wi/ojiq) = 0.719 and (a;2/wio) = 1-155, which 

are consistent with the result shown in Fig. 2. 

3. Numerical results and discussion 

In this section, the DNS technique^^' is employed to verify the theoretical results for the 
transition frequencies. The governing equations are the compressible Navier-Stokes equations 

dp ^ ^ 
— + u • Vp = -pV ■ u, 

ot 

du Vp 1 2 \ Fst 



+ u • Vu = ^ + - 2V • (/xT) - - V(/xV • u) + , 

dt P P\ 3 J P 

1^ + u • Vp = -pC|V • u, 

where p, u, p, T, p,, Fgt, and C5 denote the density, the velocity vector, the pressure, the 
deformation tensor, the viscosity coefficient, the surface tension as a volume force, and the local 
sound speed, respectively. This system of equations is solved over the whole computational 
domain divided by grids. The materials (the bubbles and a liquid surrounding them) moving 
on the computational grids are identified by a scalar function^^ depending on a pure advection 
equation whose characteristic velocity is u. The bubbles' radii and the initial center-to-center 
distance between them are set by using the same values as used for Fig. 2, that is, Riq = 5 pm, 
R20 = 9 and D{t = 0) = 20 pm. The content inside the bubbles is assumed to be an ideal 
gas with a specific heat ratio of 1.33, an equilibrium density of 1.23 kg/m^, and a viscosity of 
1.78 X 10~^ kg/(m s). The surrounding liquid is water, whose sound speed is determined by 
Cs = ^7(p + 3172.04Po)/p with p and p being the local pressure and density, respectively. 
The other parameters are set to the same values as those used previously. The axisymmetric 
coordinate (r, z) is selected for the computational domain, and the mass centers of the bubbles 
are located on the central axis of the coordinate. The grid widths are set to be constants as 
Ar = Az = 0.25 pm, and the numbers of them in the r and the z coordinates are 100 and 
320, respectively. The sound pressure, applied as the boundary condition to the pressure, 
is assumed to be in the form of Pex = PaSinut, where the amplitude Pa is fixed to 0.3Po 
and the driving frequency is selected from the frequency range around the bubbles' partial 
natural frequencies. (For the sound amplitude assumed, nonlinear effects may not completely 
be neglected especially near the natural frequencies. However, we can, unfortunately, not 
use a lower sound pressure because of a numerical problem called "parasitic currents," ^^'^^ 
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Fig. 3. DNS results: Bubble radii [(a) (j)] and corresponding positions [(a') ~ (j')] as functions of 
time for different driving frequencies. The lower lines in (a')-(j') denote the position of the smaller 
bubble. The bubbles coalesce at the time where the number of lines becomes one, see panels (i), 
(i'), (j), and (j'). 



which rise to the surface when a small bubble or drop is in (nearly) steady state where the 
amplitudes of volume and shape oscillations are very low. Effects of nonlinearity are briefly 
discussed later.) The boundary condition for the velocity is free. 

Figure 3 displays the bubbles' (mean) radii and mass centers as functions of time for 
different and figure 4 displays the sign of the force (a), determined by observing the direction 
of the bubbles' translation, and the bubbles' pulsation amplitudes (b and c). Prom these 
figures, we know that the smaller bubble has two resonance frequencies, one at a; pa a;io/0.9 
(!=a l.lwio) and the other at a; a;io/2.2 (si 0.45a;io), though the former seems to decrease 
with time because of the repulsion of the bubbles (See Fig. 1, which reveals that the highest 
transition frequency of bubble 1, causing resonance, decreases as D increases). The former 
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Fig. 4. DNS results: (a) Sign of the force, (b) pulsation amplitude of bubble 1, and (c) that of bubble 
2, as functions of the driving frequency. The amplitudes were measured for i < 5 /is (o), and for 
f < 10 yus but until the coalescence has been observed (•). The result for ui = wio/0.85, not shown 
in Fig. 3, is presented in (a); see also Fig. 5. 



resonance frequency is higher than wio, and the latter is lower than a;2o (~ O.SSwio). These 
respective resonances are obviously due to the highest and the lowest transition frequencies 
of bubble 1, both of which correspond to the natural frequencies. The same figure shows that 
the larger bubble may have a resonance frequency, at a; uiq/2.2. 

The sign of the interaction force changed twice in the frequency region considered. In 
the region between uj = a;io/0.8 and wio/O.Q, being near to but above the higher resonance 
frequency of the smaller bubble discussed above, the attractive force turns into repulsion as 
uj decreases, and, at a; (jOiq/1.2 {k, O.SSwiq) the repulsive force turns back into attraction. 
(See also Fig. 5, which shows D{t) in the cases where the deviation in it is small.) It may be 
difficult to say, using only these numerical results, that the former reversal is not due to the 
higher natural frequency of the smaller bubble, because the reversal took place near it and 
the highest transition frequency of the larger bubble is close to it. Therefore, in the following 
we will focus our attention on the latter sign reversal, which occurred at uj between cjio and 

1^20 ■ 
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Fig. 5. Timc-dcpcndcnt distances between the mass centers of the bubbles in the cases where the 
deviation in the distance is small (for lu = ujiq/0.8, CiJio/0.85, (jJio/1.2, and wio/l.S). Note that 
the result for w = wio/O.SS, not shown in Fig. 3, is added. [The low-amplitude, high-frequency 
oscillations observed in these curves may be due to a numerical error originating in the calculation 
of such a small deviation (comparable to the grid width) on a discrete computational domain.] 

The latter reversal indicates ttiat a kind of ctiaracteristic frequency stiould exist in the 
frequency region between the partial natural frequencies of the bubbles. It is evident that 
this characteristic frequency is not the resonance frequency of the larger bubble, which is, as 
already discussed, much lower. This result is in opposition to the assumption described by 
Doinikov and Zavtrak.^'^ Also, the theory for the natural frequencies (eq. (5)) cannot explain 
this reversal, because it predicts no natural frequency in the frequency region between loiq and 
UJ20- This characteristic frequency is, arguably, the second highest transition frequency of the 
smaller bubble, as was predicted by Ida.^^ As was proved theoretically in ref. 13, resonance 
response is not observed around this characteristic frequency. In order to confirm that this 
characteristic frequency is not of the larger bubble, we display in Fig. 6 the i?2-tinie and 
Pex-time curves for the area around lo = uJio/1.2. This figure shows clearly that the pulsation 
phase of the larger bubble does not reverse in this frequency region; the bubble maintains 
its out-of-phase pulsation with the external sound {i.e., the bubble's radius is large when the 
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Fig. 6. Radius of the larger bubble (solid lines) and sound pressure as functions of time, for around 
w = uj\o/1.2. 



sound amplitude is positive), although other modes, which may result from the transient, 
appear. 

Here, we discuss how this transient and a nonlinear effect act on the quantitative nature of 
the sign distribution. In Fig. 7(a), we show the sign of the force as functions of cj, determined 
by the linear theory (the solid line) and by the DNS (the circles). The positive value denotes 
the attraction, while the negative value the repulsion. A noticeable quantitative discrepancy 
can be observed between these results, whereas they are in agreement in a qualitative sense. 
We attempt here to identify what caused this discrepancy, using a nonlinear model. Mettin 
et al.^^ proposed a nonlinear coupled oscillator model for a double-bubble system, 

(^ ^Aj^i^^f^ r2 ^ f.^^A ^Ridpi^i l^rpR2^ 

= _ l + _U,, + ______(i?,i,^), (13) 



2(7 \ fRjoY'' 2a A/iRj 



with 

Pkj = fPo + -^l (=^1 -^- = -Po-Pex fori = lor2, 

based on the Keller-Miksis model^^ taking into account the viscosity and compressibility of 
the surrounding liquid with first-order accuracy. In this system of equations, the last terms of 
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(b) Nonlinear transient 
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(c) Nonlinear transient for low driving 
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Fig. 7. Comparison of the DNS results with (a) the linear harmonic solution, (b) the nonlinear 
numerical results for Pa = O.SPqi and (c) those for Pa = O.OOIPq- The circles denote the DNS 
results, where the positive value denotes the attraction while the negative value denotes the 
repulsion. 



eqs. (12) and (13) represent the radiative interaction between the bubbles. Using this model, 
Rj and Rj are calculated by the fourth-order Runge-Kutta method, and subsequently the time 
average in eq. (10) is performed to determine the sign of the force. Although the quantitative 
accuracy of this model may not be guaranteed for a small D,^^ a rough estimation of the 
influences of the transient and nonlinearity might be achieved. The physical parameters are 
the same as used for Fig. 1, and D is fixed to 20 /xm, that is, the translational motion is 
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Fig. 8. Bubble surfaces for w = uj\o/1.2 at selected times, given by the DNS technique. Similar figures 
for w = wio and w = wio/l.S can be found in our previous paper. 



neglected. The solid and the dashed lines displayed in Fig. 7(b) show 



G{T) = sgn 



.JO 



V^V2dt T 



for T = 5 /xs and 10 jis, respectively, given using the nonlinear model, where sgn[/] = 1 for 
/ > and sgn[/] = — 1 otherwise. These numerical results are in good agreement with the DNS 
results (the circles). Figure 7(c) shows the numerical results given using the nonlinear model 
with a very low driving pressure {Pa = O.OOIPq)- These results are also in good agreement 
with the DNS results, proving that, in the present case, the nonlinearity in the pulsations is 
not dominant for the sign reversals. Furthermore, it can be proved that the bubbles' shape 
oscillation is also not dominant for the sign reversal. Figure 8 shows the bubble surfaces for 
uj = LOio/1.2 at selected times. Only a small deformation of the bubbles can be found in 
this figure. (The numerical results provided previously in ref. 18, given using the same DNS 
technique, reveal that the bubbles' sphericity is well maintained even for lo = coio, whereas a 
noticeable deformation is observed for uj 0^20 ■) These results allow us to consider that the 
transient is one of the dominant origins of the noticeable quantitative discrepancy found in 
Fig. 7(a). 
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There arc some other physical factors that may also be able to cause the quantitative 
discrepancy but are not taken into account in both the linear and nonlinear coupled oscillator 
models. It is known, for example, that the translational motion of bubbles can alter the bub- 
bles' pulsation. As has been proved by several researchers,^^' ^"^'^^ if the translational motion 
is taken into consideration in deriving a theoretical model, high-order terms appear in the 
equations of radial motion, in which not only the bubbles' radii but also their translational 
velocities are involved. Also, in certain cases the viscosity of the surrounding liquid can alter 
the magnitude and sign of the interaction force, because acoustic streaming is induced around 
the bubbles. Investigating the influences of those factors on the quantitative discrepancy 
using a higher-order model would be an interesting and important subject. 

4. Concluding remarks 

In summary, we have verified the recent theoretical results regarding the transition frequen- 
cies of two acoustically interacting bubbles^^ and the sign reversal of the secondary Bjerknes 
force,^^ using a DNS technique. ^^'^^ The present numerical results, given by the DNS tech- 
nique, support those theoretical results at least in a qualitative sense. The most important 
point validated by DNS is that at least the sign reversal occurring when the driving frequency 
stays between the bubbles' partial natural frequencies is obviously not due to the natural (or 
the resonance) frequencies of the double-bubble system. This conclusion is in opposition to the 
previous explanation described by Doinikov and Zavtrak,^'^ but is consistent with the most 
recent interpretation by Ida^^ described based on analyses of the transition frequencies,^^' 
thus validating the assertion that the transition frequencies not corresponding to the natural 
frequencies exist and that the notion "transition frequency" is useful for understanding the 
sign reversal of the force. 

Acknowledgment 

This work was supported by the Ministry of Education, Culture, Sports, Science, and 
Technology of Japan (Monbu-Kagaku-Sho) under an IT research program "Frontier Simula- 
tion Software for Industrial Science." 



16/18 



J. Phys. Soc. Jpn. Full Paper 

References 

L. A. Crum: J. Acoust. Soc. Am. 57 (1975) 1363. 
A. Prosperetti: Ultrasonics 22 (1984) 115. 

W. Lauterborn, T. Kurz, R. Mettin and C. D. Ohl: Adv. Chem. Phys. 110 (1999) 295. 
E. A. Zabolotskaya: Sov. Phys. Acoust. 30 (1984) 365. 
A. A. Doinikov and S. T. Zavtrak: Phys. Fluids 7 (1995) 1923. 
A. A. Doinikov and S. T. Zavtrak: J. Acoust. Soc. Am. 99 (1996) 3849. 
Y. A. Kobelev, L. A. Ostrovskii and A. M. Sutin: JETP Lett. 30 (1979) 395. 
P. L. Marston, E. H. Trinh, J. Depew and T. J. Asaki, in: Bubble Dynamics and Interface Phe- 
nomena, edited by J. R. Blake et al. (Kluwer Academic, Dordrecht, 1994), pp. 343-353. 
P. C. Duineveld: J. Acoust. Soc. Am. 99 (1996) 622. 

P. A. Dayton, K. E. Morgan, A. L. Klibanov, G. Brandenburgcr, K. R. Nightingale and K. W. 
Ferrara: IEEE Trans. Ultrason. Ferroelect. & Freq. Control 44 (1997) 1264. 
P. Dayton, A. KUbanov, G. Brandenburger and K. Ferrara: Ultrasound Med. Biol. 25 (1999) 1195. 
A. Shima: Trans. ASME J. Basic Eng. 93 (1971) 426. 
M. Ida: Phys. Lett. A 297 (2002) 210. 
M. Ida: Phys. Rev. E 67 (2003) 056617. 
M. Ida: J. Phys. Soc. Jpn. 71 (2002) 1214. 

T. J. Matula, S. M. Cordry, R. A. Roy and L. A. Crum: J. Acoust. Soc. Am. 102 (1997) 1522. 
1. Akhatov, R. Mettin, C. D. Ohl, U. Parlitz and W. Lauterborn: Phys. Rev. E 55 (1997) 3747. 
M. Ida: Comput. Phys. Commun. 150 (2003) 300; Erratum in 150 (2003) 323. 
M. Ida: Comput. Phys. Commun. 132 (2000) 44. 
T. Yabe and P. Y. Wang: J. Phys. Soc. Jpn. 60 (1991) 2105. 
S. Ito: 43rd Nat. Cong, of Theor. & Appl. Mech. (1994) p.311 [in Japanese]. 
M. Ida and Y. Yamakoshi: Jpn. J. Appl. Phys. 40 (2001) 3846. 
J. U. Brackbill, D. B. Kothe and C. Zemach: J. Comput. Phys. 100 (1992) 335. 
T. G. Leighton: The Acoustic Bubble (Academic Press, London, 1994), p.291. 
C. Devin: J. Acoust. Soc. Am. 31 (1959) 1654. 
A. Prosperetti: Ultrasonics 22 (1984) 69. 

S. Hilgenfeldt, D. Lohse and M. Zomack: Eur. Phy. J. B 4 (1998) 247. 
C. Feuillade: J. Acoust. Soc. Am. 98 (1995) 1178. 
Z. Ye and C. Feuillade: J. Acoust. Soc. Am. 102 (1997) 798. 

A. Harkin, T. J. Kaper and A. Nadim: J. Fluid Mech. 445 (2001) 377. 
H. Takahira, S. Fujikawa and T. Akamatsu: JSME Int. J. Ser. II 32 (1989) 163. 
P.-Y. Hsiao, M. Devaud and J.-C. Bacri: Eur. Phys. J. E 4 (2001) 5. 

B. Lafaurie, C. Nardone, R. Scardovelli, S. Zaleski and G. Zanetti: J. Comput. Phys. 113 (1994) 
134. 

S. Popinet and S. Zaleski: Int. J. Numer. Methods Fluids 30 (1999) 775. 

R. Mettin, I. Akhatov, U. Parlitz, C. D. Ohl and W. Lauterborn: Phys. Rev. E 56 (1997) 2924. 
J. B. Keller and M. Miksis: J. Acoust. Soc. Am. 68 (1980) 628. 
H. Oguz and A. Prosperetti: J. Fluid Mech. 218 (1990) 143. 
H. Takahira, T. Akamatsu and S. Fujikawa: JSME Int. J. Ser. B 37 (1994) 297. 



17/18 



J. Phys. Soc. Jpn. Full Paper 

39) A. A. Doinikov: J. Acoust. Soc. Am. 106 (1999) 3305. 



18/18 



